Pairing mechanism in Fe pnictide superconductors 
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By applying an exact unitary transformation to a two-band hamiltonian which also includes the 
effects due to large pnictogen polarizabilities, we show that an attractive spin-mediated Hubbard 
term appears in the d xz , d yz nearest-neighbour channel. This pairing mechanism implies a singlet 
superconducting order parameter in iron pnictides. 

PACS numbers: 71.10.-w, 74.20.-z, 74.20.Mn 



The recent observation of high temperature supercon- 
ductivity in Fe pnictide compounds [l[ has generated a 
huge interest and a large research effort. As a result four 
families of compounds have emerged having similar prop- 
erties and exhibiting superconductivity. Namely, chalco- 
genides (FeSe and FeTe) and three families of pnictides 
(represented by LaOFeAs, BaFe2As2, and LiFeAs). 

Band structure calculations have revealed that (i) 
electron-phonon coupling alone cannot account for the 
high values of T c Q , (ii) nesting between the hole Fermi 
surfaces at zone center and the electron Fermi surfaces 
around the zone corner may play a significant role both 
in magnetism and superconductivity [3]. These results 
have prompted the idea that the pairing glue is provided 
by spin fluctuations exchange between electrons in dif- 
ferent bands [H, However, NMR studies measuring 
the temperature dependence of nuclear spin-lattice re- 
laxation rate (1/Ti) suggest that spin fluctuations may 
not be important in iron-pnictide superconductivity [5|. 
Thus other pairing mechanisms have to be responsible 
for the unusually high superconducting critical transition 
temperature. 

In the following we will present the results of an ex- 
act unitary (canonical) transformation (UT) by which it 
can be shown rigorously that an attractive spin-mediated 
Hubbard term appears in the d xz , d yz nearest- neighbour 
channel. The mechanism is similar to a spin double- 
exchange type and hence has origins of a kinetic mecha- 
nism. 

We start with the by now familiar two-band model 
description of the FeAs superconductors Q. We take a 
two-dimensional square lattice with d xz and d yz orbitals 
per site. For simplicity we label the d xz and d yz orbitals 
by a and b, respectively. The kinetic energy component 
can be expressed as: 



Hkin — 



E M a tA* + h - c -) 



(1) 



ij,cr a,f3—a,b 



Here, the site indices, i and j, run over nearest and next- 
nearest neighbors. The effective hoppings, t a p, have con- 
tributions from both direct Fe-Fe and Fe-As-Fe processes. 



If the simplified description of a Slater-Koster formalism 
is adopted, then t a b, the matrix element between d xz 
and d yz orbitals on nearest neighbor sites, is zero because 
the formalism assumes perfect spherical symmetry. How- 
ever, in models that use all five Fe 3d orbitals, t a b is non- 
zero 0, H|. In fact, according to Ref. @], t ab = 0.54 eV 
is the largest Fe-Fe hopping integral. Consequently, we 
take the view that the two-band model can be extracted 
from a more general model by focusing on its relevant 
2x2 block. The parameter t a b is crucial because, as will 
be discussed below, it will be at the basis of the UT (see, 
Eq. ((4])). The other hopping parameters enter into our 
UT formulation through the bandwidth, e, which we take 
as 2 eV (KEHH. 

-Hint = 5Zi H\ contains only on-site contributions, with 



Hi = J2 a U a ni, a ,im, a ,i + (U 1 - J/2) J2i ni, a ni,b 
-2JS la • Si, 6 + J'(ai,T°U 6 i -l b i d + h - c -) ' 



(2) 



As in Eq. ([I]), a = a,b labels d xz ,d yz orbitals. Si, a 
ini, a ,cr) is the spin (density) in orbital a at site i. Follow- 
ing Refs. Q, we used U' = U — 2 J and the pair hopping 
term strength J' = J, where J is the Hund coupling. 

The onsite Hubbard terms are obviously equal, U a = 
Ub = U, with U chosen between 3.0 - 4.5 eV. Although in 
the literature there is no agreement on whether the iron 
pnictides should be considered as strongly correlated or, 
at most, moderately correlated materials, there is more 
consensus regarding the strength of U to be 4 eV fioj ]. 
3.5 eV [IH or even smaller U > 2 eV 0. In our UT 
approach described below, U' and J give rise to no new 
physics, although they do render the transformation alge- 
bra very complicated. Thus, only at the end their effect 
on the final result will be briefly discussed. 

The third and last contribution that we include in 
our starting hamiltonian is the polarizability effect. As 
pointed out already in Ref. [3] the polarizability effects 
in iron pnictides are much larger than in copper-base 
high T c 's. In fact, it is known [14| that the ions As 3 ~ 
and Se 2 ~ have huge polarizabilities, mainly due to the 
fact that their volume is very large. Hence, whenever an 
iron site is charged due to electrons hopping to or from it, 



the surrounding As or Se atoms will easily get polarized, 
an effect we need to capture. 

The effect of As (or Se) polarizability on iron sites can 
be described by writing [1 51 ] the hamiltonian term first 
introduced in Ref. [l6| per iron site as g n (jfts + s^p). 
Here the possible excitations of an As electron from 4p 
to unoccupied 5s are taken into account with an effective 
coupling, denoted by g, due to a charge on Fe. The n — 
rij + ni notation (n — n a or n b ) is identical to the one 
used in Eq. |(2~J) and describes Fe 3c? electrons. 

The angle dependence between the As p and s or- 
bitals was extensively studied already [17[, hence we 
only consider a simplified version of the hamiltonian 
taking into account only the mean total polarizability 
per bond (lEL [l6| by a hamiltonian term of the form: 
H poi = p J2(ij)( n i ~ where P = g(j?h + s^p) will 
be the measure of As average polarizability, with P es- 
timated to be around 2eV 17j. In this way our starting 
hamiltonian is H^ n + H lni + -ff po i. 

One of the most common tools used in theoretical 
physics is perturbation theory. Here, we present a pertur- 
bation theory which we solve exactly using a UT. There 
are two reasons for using UT [l^|: i) the belief that the 
transformed Hamiltonian is "simpler" in the sense that 
it is "more diagonal" and ii) the desire to gain a deeper 
physical insight into the problem, given that the trans- 
formed Hamiltonian may reveal the appropriate indepen- 
dent subsystems. Our scope is to pursue the latter. 

For the UT to work, we separate out from -ffkin the di- 
agonal hopping term, H ab = t ab Y,{ij), a { c H.,J ) i,a + n - c -)- 
Hence we can write H^ n = Hq + H abl where Hq contains 
the remaining terms other than t ab . To be more trans- 
parent, we re- write our starting hamiltonian as H + H a {>, 
where H = Hq + H mt + H po \. H will be the zeroth-order 
hamiltonian and t a b the perturbation. 

1^ . 2o| that in fourth order perturbation 



It is known 



an attractive Hubbard type carrier-carrier interaction ap- 
pears, parallel to spin-carrier, spin-spin, carrier-carrier- 



spin, etc. terms |2l|, [22j. All these terms are spin me- 



diated, of superexchange type [3, 0], and hence with 
physical origins of a kinetic mechanism [23| . The attrac- 
tive carrier-carrier terms that we are interested in appea r 
either onsite or between nearest-neighbor sites [l9|, |20j . 
However, performing a perturbation up to fourth order 
is not enou gh t o draw a final conclusion, as was pointed 
out in Ref. [201 ] . Depending upon the value of input pa- 
rameters, the second, sixth, etc., order terms may well 
give repulsive interactions, while the fourth, eighth, etc., 
attractive ones 22]. Thus, higher contributions need to 
be calculated in order to verify the convergence. 

Hence, initially we check for the convergence of higher 
order terms for the particular case of Fe pnictides. In a 
standard UT [3, dS [2(1 HH, E3| the transformed Hamil- 
tonian e s (H + H a b)e~ s is identical to H + X]^li[V n ' — 



'•h o.o 




FIG. 1: The calculated superexchange interaction J a t in dif- 
ferent orders n, as defined in Eq. @. We show the total 

contributions for each order, i.e., 2~Z™=i[l/ n - — V( n + 1)'] H n , 
see also Ref. [24[. For input parameters we used e = 2eV, P 
= 2eV and U = 4.5eV. 



l/(n+l)!]-ff n (henceforth, denotes a UT result), where 



n times 



H,, 



[S,[S,[S, ■■■ [S,H ab 



(3) 



for an S which satisfies H a b + [S, H] = 0. Using the 
notations: C x = l/(s + P), C 2 = l/(e + P + U b ) - d, 
C 3 = l/(e + P -U a ) - d, and C 4 = l/(e + P + U b - 
U a ) ~ l/(e + P + U b ) - l/(e + P-U a )+Cx, the unitary 
operator S is: 

S = -UE<ij) )CT (Ci +C2nl_ a + C 3 nZ_ (7 

+C i n?_ a nl_ a ) (4^. ct - h.c.) . (4) 

The first order UT (i.e., second order perturbation the- 
ory [25I ] ) c an be readily performed. This result is well 
known [20j, [22|, and consequently not pursued further. 
In the second and all other 'even' order transformations, 
H a b is recovered structurally with the addition of all 
possible correlated hopping terms [11, H(|: H n=evcn = 



Z)(ij),tr[*o6+*o. 



b'H,- 



ab' l i.- 



»'■' 



Consequently, the weight of the original t ab is redis- 
tributed in every order of the transformation among t ab , 
t a ab , t^ b and t" b . Because of the smallness of U in Fe pnic- 
tides, as the UT is performed to higher orders weight is 
shifted in and out of the correlated hopping terms gen- 
crating an oscillatory behavior. These oscillations are 
well-known 3, 2^, 23] m standard perturbation theory. 
As an example, the Heisenberg superexchange term J ab 
is shown in Fig. [1] for high orders. Hence, the problem 
with applying standard UT is that unless U is very large, 
the perturbation series will not converge. 

To overcome this problem, a completely new approach 
is needed to handle Fe pnictide case, i.e., to perform an 
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FIG. 2: (Color online) The calculated nearest neighbor, U a b 
(top), and the total onsite, U (initial value) + U (calculated 
value) (bottom), Hubbard coefficients as a function of tab- 
Input parameters are s = 2eV, U = 3eV [part (a)] and U — 
4.5eV [part (b)]. Continuous (black) curves are P = 1.5eV, 
while long-dashed (red) curves are for P — 3eV. 



exact UT. We have chosen to perform such a transforma- 
tion by eliminating three consecutive even order terms 
simultaneously (the n — 0, n — 2 and n — 4 order terms) 
from Eq. ([3]) via a new unitary operator S such that 



H ab +[S,H] + -[S, [S,H ab ]} 







(5) 



This guarantees that the transformation cannot be con- 
tinued to higher orders, since all a <-> b hopping pro- 
cesses are eliminated, as t ab and the correlated hopping 
terms, t a ab , t a [ b , tf* b are strictly zero. Accordingly, the 
transformed hamiltonian 



e s (H + H ab )e- 



H + [<S, H ab 



[S,[S,H], (6) 



is exact in a strict mathematical sense [18|. Eq. ^ is 
merely a self-consistent mixing 18J of the n=l and n=3 
order standard UT, i.e., of a standard 2nd and 4th order 



perturbation theory 25|. This guarantees that we can 
capture the whole spectrum of the carrier-carrier inter- 
action terms • 

The solution of Eq. has the same algebraic form 
as Eq. (0} with new coefficients C%, C 2 , C3 and C4 as 
unknown parameters. These are determined from Eq. 
(|5|), which gives the following systems of equations: 
^ = U[cos(2UCi) + (e + P)d sin{2t ab C 1 )/(2t ab C 1 )} 
C = Kb^iz + P)( M 2 - v 2 ) - U a p 2 - U b v 2 \^2r) + 
8^ 6 m0* " ^(2t) + ttl b nv 2 (e + P)c^{t) + 2t\ b v{p. + 
v)^{t) - t ab cos(2UCi) - (l/2)(e + P) sm(2i ah Ci), where 
H = & + C 2 , v = d + C 3 , t = t ab [2{p 2 + v 2 )] 1 / 2 and the 
notations (fr(x) = sin(x)/x 3 , ip(x) — cos(x)/x 3 were used. 
For t^ b we have to exchange U a with U b . Finally, the last 
of the correlated hopping terms is: 

ttl = -8^0" ~ v)[(e + P)(m 2 - v 2 ) ~ C/aM 2 - 
Ukv 2 }^) - m 3 ab (» - v) 2 ^(2t) - 2t 3 abf iv(fi + v) \{e + P) - 



U a + U b ]^{r) - 2tl b {n+v) 2 i>{T) +t ab {Cl + p 2 ) cos(2t abP ) + 
±(e + P)(l + Cl/p 3 ) sin(2Up) + U(l - C 2 /p 2 )[l + (e + 
P)Ci], where p = C± + C 2 + C 3 + C4 and note that these 
equations can also be used if U a ^ U b . 

We calculated numerically the parameters C\, C2, C3 
and Ci from the above system of equations by setting t ab , 
t® b , l^ ab and t a a h b to zero. The new UT, as any standard per- 
turbation theory, gives several spin-spin, carrier-carrier, 
spin-carrier onsite, nearest neighbor, next-nearest neigh- 
bor, etc. terms. From all these we only calculated ex- 
plicitly the onsite, U^2 i ni^^ni a^ (a = a, 6) and near- 
est neighbor U ab Y^,n j\ ™i.a,T n j,b,i Hubbard interactions. 
The results are plotted in Fig. H as a function of t ab , 
and in Fig. [3] as a function of the polarizability, P. 

As can be seen in Fig. ® a nearest-neighbor inter- 
orbital attractive Hubbard interaction U ab appears due 
to a spin-mediated superexchange type mechanism 
I20II . Being of kinetic origin, the spins are not excited 
[23j, only virtual excitations of onsite singlets occur. The 
situation is similar to the phonon-mediated attraction: at 
temperatures well below the Debye scale real phonons are 
never excited, yet they provide an attraction mechanism 
between electrons. 

The values of U ab in Fig. [2] can be well approxi- 
mated for small t ab with a line: U ab ~ — const t abl 
with the constant being 4£(C+ - C-) {sin(2[2£ 2 (C 2 + 
£)] 1/2 )/(2[2£ 2 (# +C 2 )] 1/2 ) + [cos(2[2£ 2 (a + C 2 )] 1/2 )- 
1]/(2[2£ 2 (C 2 + C 2 )] 1/2 ) 2 }, where C± = l/[e + P ± U), 
£ = ?7|(e + P) 2 -t/ 2 |/[(e + P) 2 + C/ 2 ] 1 / 2 and the numerical 
coefficient rj is 0.392 for e + P < U, or 1.178 otherwise. 

The new hamiltonian term, Hp Q \, which gives a mea- 
sure of the polarizability is not crucial in obtaining this 
attractive interaction. Any finite value of bandwidth 
alone suffices for attraction. In fact, the UT expressions 
depend only on the sum e + P, and not separately on 
either. Thus, polarizability effectively acts as a band- 
width. So, Fig. [3] covers a wider range for P to allow 
for cases with different bandwidth and P values. From 
Fig. [3l upper panel we observe that the attractive U ab 
is enhanced rather strongly by P for e + P ^ U. The 
strongest attraction occurs around the values of P satis- 
fying e + P ~ U. Note that, this regime is not accessible 
by standard perturbation theory, but it works well in an 
exact UT approach. 

An additional effect due to Hp Q \ is the modified screen- 
ing of the on-site Coulomb repulsion U in the presence 
of P. In the lower panel of Fig. [3] the net repulsion, 
U + U, is seen to deviate significantly from a roughly 
constant value (< 3 eV for U — 3 eV, « 4 eV for U = 4.5 
eV) only in an interval given by \e + P — U\ < 0.5 eV 
= t ab . Normally, ligand polarizabilities are expected to 
increase the screening (smaller net U values), however, 
when the double exchange is the strongest, it restricts 
the hopping and as such screening is reduced, similarly 
to, e.g., manganites [14|. U' and J from Eq. have 
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FIG. 3: (Color online) The same Hubbard coefficients as in 
Fig. [2] as a function of the polarizability, P, for e = 2eV and 
tab = 0.5eV. U = 3eV, continuous (black) curves and 4.5eV 
long dashed (red) curves, respectively. For the lower curves, 
the left (right) axis is for U = 3eV (4.5eV). 



been neglected up to now. However, their effects in our 
UT can be accounted for by substituting everywhere U 
with U + J(l - U'/U)/2. Using the known U' and J 
values [1, d, 0, [13, El, E2] this only causes a less than 
5% change in U and hence their effect is minor. 

In conclusion we have shown that the Fe pnictide fam- 
ily of alloys exhibit a rather strong pairing interaction, 
due to a spin-mediated, superexchange type mechanism. 
There is a quite strong enhancement of the attraction 
due to the polarizability. However we do not observe 
polaron and bipolaron formation as in Ref. [171 ]. proba- 
bly because we are using realistic onsite Hubbard values, 
instead of U > 8 eV as used in Ref. [13] ■ Our pair- 
ing mechanism will give rise to a singlet superconducting 
order parameter. However, in a two-band system the 
self-consistent solution of the gap equations always has 
a symmetric (recently denoted as s++) and asymmetric 
(s-| ) solution 26]. In the parameter regime of supercon- 
ducting Cr alloys, at different concentration of electrons 
and holes, the asymmetric solution wins [27|. In clos- 
ing we remark that non-phononic mechanisms support 
[28j negative isotope effect [§9|. The standard isotope 
effect measured [3fJ| for both superconductivity and itin- 
erant antiferromagnetism in Fe pnictides would require 
phonon input, which may contradict an electronic mech- 
anism. However, this may not necessarily be the case, as 
phonons can act as random impurity potentials in certain 
cases for systems exhibiting electron and hole Fermi sur- 
faces [3l| . The pair-breaking effect originating from this 
also explains why Cr has an isotope effect in TV. This 
actually can be very large; the phonons can decrease TV 
by as much as 70% from its value in their absence [32 ]. 
Another interesting consequence of the pair-breaking by 
virtual or thermal phonons [33] is that the ratio 2A/ksT c 
can be much larger than the BCS value of 3.53. 



MG thanks L. Hozoi and P. Fulde for fruitful discus- 
sions. He also acknowledges the financial support of the 
Max Planck Institute for the Physics of Complex Sys- 
tems, Dresden, where most of this work was carried out. 



[1] Y. Kamihara, et al., J. Am. Chem. Soc. 130, 3296 (2008). 
[2] L. Boeri, et al, Phys. Rev. Lett. 101, 026403 (2008). 
[3] I. I. Mazin, et al., Phys. Rev. Lett. 101, 057003 (2008); 

A. V. Chubukov, et ah, Phys. Rev. B78, 134512 (2008). 
[4] K. Kuroki, et al, Phys. Rev. Lett. 101, 087004 (2008). 
[5] N. Terasaki, et al., J. Phys. Soc. Jpn. 78, 013701 (2009); 
K. Ishida, et al., J. Phys. Soc. Jpn. 78, 062001 (2009); 
Y. Nakai, et ah, Phys. Rev. B79, 212506 (2009). 
[6] S. Raghu, et al., Phys. Rev. B77, 220503(R) (2008); M. 
Daghofer, et al., Phys. Rev. Lett. 101, 237004 (2008); A. 
Moreo, et al, Phys. Rev. B79, 134502 (2009). 
[7] J.C. Slater and G.F. Koster, Phys. Rev. 94, 1498 (1954). 
[8] E. Manousakis, et al., Phys. Rev. B78, 205112 (2008). 
[9] C. Cao, et al, Phys. Rev. B77, 220506(R) (2008); K. 
Nakamura, et ah, J. Phys. Soc. Jpn. 77, 093711 (2008). 
[10] K. Haule, et al, Phys. Rev. Lett. 100, 226402 (2008). 
[11] V.I. Anisimov, et ah, Physica C469, 442 (2009). 
[12] W.L. Yang et ah, Phys. Rev. B 80, 014508 (2009). 
[13] H. Takahashi, et al. Nature 453, 376 (2008). 
[14] P. Fulde, Electron Correlations in Molecules and Solids, 

Springer- Verlag, 1991. 
[15] M. B. J. Meinders, et al, Phys. Rev. B52, 2484 (1995). 
[16] D. K. G. de Boer, et al, Phys. Rev. B29, 4401 (1984). 
[17] G. A. Sawatzky, et ah, Europhys. Lett. 86, 17006 (2009); 

M. Berciu, et al., Phys. Rev. B79, 214507 (2009). 
[18] M. Wagner, Unitary Transformations in Solid State 

Physics, North-Holland, 1986. 
[19] V. J. Emery, Phys. Rev. Lett. 58, 2794 (1987); V. J. 

Emery and G. Reiter, Phys. Rev. B38, 4547 (1988). 
[20] J. Zaanen and A. M. Oles, Phys. Rev. B 37, 9423 (1988). 
[21] A. H. MacDonald, et al, Phys. Rev. B37, 9753 (1988); E. 
Kolley and W. Kolley, Phys. Stat. Sol. 157, 399 (1990); 
E. Kolley, et al, J. Phys.: Condens. Matter 4, 3517 
(1992); H. Eskes and J. H. Jefferson, Phys. Rev. B48, 
9788 (1993). 

[22] M. Gulacsi and R. Chan, Jour. Supercond. 14, 651 
(2001); E. Muller-Hartmann and A. Reischl, Eur. Phys. 
J. B28, 173 (2002); R. Chan and M. Gulacsi, Phil. Mag. 
84, 1265 (2004). 

[23] P. W. Anderson, Phys. Rev. 115, 2 (1959). 

[24] As the number of terms (2 1 +™+(™- 3 )/ 2 5(™- 1)/2 , n > 3) 
which contribute to J a t, in each individual odd order can 
reach tens of billions, we used a purpose built algebraic 
manipulation program to calculate higher orders. 

[25] The nth order UT transformation of the Eq. © type 
corresponds to an (n + l)th order standard perturbation 
theory [l8| . 

[26] A. I. Rusianov, et al, Sov. Phys. - JETP 38, 991 (1974). 
[27] M. Gulacsi and Zs. Gulacsi, Phys. Rev. B33, 6147 (1986). 
[28] D.R. Penn, et al, Ph ys. Rev. B47, 5528 (1993). 
[29] P. M. Shirage, et al., arXi v: 0903. 35151 (2009). 
[30] R. H. Liu, et al, Nature 459, 64 (2009). 
[31] A. S. Barker, Jr., et ah, Phys. Rev. Lett. 20, 384 (1968). 
[32] E. W. Fenton and C. R. Leavens, J. Phys. F: Metal Phys. 
10, 1853 (1980). 



